Europalsches Patentamt 
European Patent Office 
Office europeen des brevets 



© Publication numder 



0 368 436 

A2 



© 



EUROPEAN PATENT APPLICATION 



© Application number 89304479.2 
© Date of Wing: 04.05.89 



© int. ci * G01N 21/47, G01N 2i/01, 
A61B 5/00 



© Priority: 08.11.88 GB 8826104 

© Oate of publication of application: 
16.05.90 Bulletin 90720 

® Designated Contracting States: 
DE GB 



© Applicant HAMAWATSU PHOTONICS K.K. 
1126*1 Ichlno-csso Hamamatsu-shl 
Sh izuoka-ken<JF) 

© Inventor: Van der Zee, Plet of Univ. ColL 
London Dept. Med. 

Physics and Bfc-Eng. 1st Floor Shropshire 
House 

n-20 Capper Sr«et London WC1E SJA(OB) 
Inventor: Delpy. Qavld T. of Univ. ColL London 
Dept Med. 

Physics and B9c*€ng. 1st Floor Shropshire 
House 

11-20 Capper Sceet London WC"!E 6JA(GB) 
Inventor Arrfdg*. Simon of Univ. ColL London 
Dept Med. 

Physics and Bfc-Eng. 1st Floor Shropshire 
House 

11-20 Capper Swet London WC1E 6JA(GB) 



© Representative: Beckham. Stephen Nell et al 
GILL JENNINGS & EVERY 53-64 Chancery 
Lane 

London WC2A THN(GB) 



0 Phantom having similar optical characteristics to living tissues. 

© A phantom having similar light scattering characteristics to Irving issues of interest is obtained by 
suspending in a selected medium (2) a mixture of particles (3) of differens-*Sued radius r k at respective values 
C»J of concentration C fc which are determined on the basis of the following matrw equation: 
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UJ where P^fl.) is a scattering Intensity of light with the living tissues at a scattering angle 9, and P(MiJ is a 
scattering intensity with the particles having the radius r„ at the scattenrc; angle «,. 
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The present invention rotates to a phantom having optical characteristics, in particular light-scattering 
characteristics, that are similar to those of living tissues. 

Cerebral hemorrhage and inadequate supply of oxygen to brain cells have been either the primary 
cause of the death of newborn babies in an ICU (Intensive Care Unit) or one of the causes that make them 
handicapped throughout their life. Bleeding in the brain can be readily discovered by ultrasonic diagnosis or 
by other imaging techniques. However, hypoxia of brain cells, as with cysts in porencephalic patients and 
other deficiencies of brain tissues, cannot be detected until after about 3 weeks have passed since the 
occurrence of the damage. Hypoxia in the brain cells of newborn babies has even become a social concern 
for medical institutions in developed countries. In Great Britain, intranatal deaths resulting from asphyxia 
account for one third of the causes of natal death. In other words, three out of 1 .000 newborn babies die of 
suffocation and at least an equal number of babies who survive will be handicapped on account of asphyxia. 

Under these medical and social circumstances, a need has arisen for the development of a diagnostic 
apparatus that allows direct measurement of oxygen supply to brain cells. Tissues and bones of human 
bodies exhibit a fairly good transparency to light in the near infrared region. Furthermore, in this wavelength 
range, hemoglobin and cytochrome which ere vectors of oxygen transport through the body undergo a 
change in their absorption spectra depending upon whether they are oxidized or reduced. Therefore, the 
change in the amount of oxygen supply to the brain can be detected by irradiating the head with near 
infrared light in the wavelength range of 700 - 900 nm and measuring the change in the spectrum of light 
that has been transmitted through the head. 

Various diagnostic devices have so far been proposed on the basis of the principle described above 
and their operation consists essentially of guiding light from a near infrared laser diode to the head via fiber 
optics and detecting tht i emitted light from the head (see. for example. U.S. Patent Nos. 4.223.680 and 
4.281 .645). These devices are adapted to measure the overall (average) change in the quantity of oxygen in 
the brain. Techniques for determining the distribution of oxygen in the brain have also been investigated 
and some of them are discussed in "Bioinstrumentation Using Light - Roads to Optical CT* in O plus E. 
May 1987 to April 1988. 

One of the prerequisites for the development of devices of the type described above is the use of a 
phantom, or a volume of material that has optical characteristics (absorption and scattering) equivalent to 
those of brain tissues, and this phantom is required to be stable and to ensure reproducible results. 

An oiltn-water emulsion has been conventionally used as a phantom but this is not highly stable and 
suffers the disadvantage of experiencing a change ir. particle size with time. A further problem is that the 
particles in the emulsion often become too small to realize optical characteristics similar to those of living 
tissues of interest 

According to the present invention, a phantom has simitar optical characteristics to living tissues with 
respect to light-scattering characteristics such as the scattering distribution and has a structure that in a 
selected medium is suspended a mixture of particles of different-valued radius r k (k is a positive integer of 1 
to m) at respective values of concentration C* which are determined on the basis of the following matrix 
equation: 
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where P««(9.) is a scattering intensity with the living tissues at a scattering angle of * (i is a positive 
integer of 1 to n) and P (fl^d) is a scattering intensity with a particle having the radius M. 

According to the present invention, a phantom having similar optical characteristics to living tissues of 
interest can be produced by merely incorporating different-sized particles in a medium. Particularly high 
stability and reproducibility can be reafised by selecting particles whose specific gravity is generally the 
same as that of the medium (this may be achieved by mixing polystyrene particles with water for example). 

Examples of phantoms in accordance with this invention will now be described with reference to the 



2 



BP 0 368 436 A2 



JO 



15 



accompanying drawings in which:- 

Rg. 1 shows how to determine a distribution of fight scattering on the basis of measured daa: 
Fig. 2 shows specific procedures of forming a phantom; 

Fig. 3 is a graph showing the scattering characteristics of a phantom composed of two kinds of 
polystyrene particles; and. 

Fig. 4 and S are diagrams i'tustrating two typical uses of t?ie phantom of the present invention. 
A phantom of the present invention is prepared by suspending different-sized particles in a medium. 
The particles to be suspended desirably have a specific gravity that is generally the same as that of the 
medium. This condition may be satisfied by suspending polystyrene particles in water. The wavelength 
range in which optical characteristics of interest are to be measured is desirably that of near infrared 
radiation (e.g. 700 - 900 nm) and the absorption of light by polystyrene particles is negligible in this 
wavelength range. 

The phantom of the present invention is formed by suspending a mixture of particles, say, pcCystyrene 
particles, of different-valued radius r„ (k is a positive integer of 1 to m) at respective values of concentration 
C k in a selected medium, say. water. The respective values of C k are determined as follows, if the 
measured scattering distribution with living tissues at a scattering angle of 6. (i is a positive integer of 1 to n) 
is written ss P«»„(S ( ) and if the calculated or measured scattering distribution with polystyrene particles 
having the racSus r k is written as P(9 t .r0. then the concentration. Ck. of polystyrene particles having a 
specific value of radius r„ can be determined on the basis of the following matrix equation: 
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The scattering distribution. Pm««(9i). given on the left side of equation (1) can be obtained by measuring 
& the scattering characteristics of the tissues of a living organ such as the brain illuminated with incident light 
If a profile as shown in Fig. 1 is obtained, which represents the dependenco of the intensity otf scattered 
light on the scattering angle, the intensities of scattered light at scattering angles e i . $ it and B m correspond 
to P^Bi). P-^,(fli) and P«*«(M. respectively. 

The scataring distribution, P(«Vh). given in the first term of the right side of equation <i> «s that of 
» polystyrene particles having a radius of r*. In general, P(9.M is calculated by known Mie theory of light 
scattering but it may also be determined by actual measurements, in this latter case, the scattering 
characteristics of polystyrene particles are measured for each value of the radius r h and Pts.rJ is 
determined from the obtained data. 

On the of these measured and calculated values, operations are performed on the matrix equation , 
40 (1 ) to determne the concentration of polystyrene particles. Ck. for each value of tr* radius r». Operations 
with the matnat equation (1) may be performed by the SVD (single value decomposition) method. 

Specific procedures of forming the phantom described above are explained hereinafter with reference to 
Fig. 2. 

First. living tissues of the brain are prepared as a sample and illuminated with near infrared fight (783 
45 nm ) to investigate the characteristics of light scattering in the sample. The results are as shown graphically 
in Rg. 2(a). with the scattering angle (ft) being plotted in degrees on the horizontal axis and the density of 
scattered Bght being plotted on a logarithmic scale on the vertical axis. This complete* the step of 
determining experimentally the scattering distribution P^^,(Si). 

In the next step, the scattering distribution P{9ij^ is calculated by We theory of tight sc atterrxj On the 
so basis of both the results of this calculation and the measured values shown In Fig. 2(a). the concentration o? 
polystyrene particles. C k . is calculated for each value of the radius r fc by equation (1). As for the scattenng 
angle this calculation is performed over the range of 0-1 80 " with 1 " increments, and as for the radius r». 
calculation is carried out over the range of 0 - 0.8 urn with 0.04 am increments. The results are as shewn 
graphically in Fig. 2(b). with the horizontal axis plotting the radius of polystyrene particl es. r» . a nd the 
93 vertical axis plotting relative concentration (on a linear scale). One can readily see the presence of two 
peaks in the graph as denoted by A and B. 

The scaaering distribution. P(e lt r^. as determined by Mie theory of Bght scattering is shown graphically 
in Fig. 2(c) w*h the scattering angle plotted on the horizontal axis. The horizontal axis of the graph shown in 
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Fo 2fc) corresponds to the scattering angle and has been normalized ;n accordance wd the SVD method. 

* When polystyrene particles were suspended in water at respective values of the concentration C* 
destined for each value of the radius r k (see Fig. 2(b)). a phantom was obtained that had tight-scattering 
ci%a-sc ten sties simitar to those of the living tissues. In other words, when the intensity of scattered light was 
s cacuiated on the basis of the data of concentration C* shown in Fig. 2(b) for each value of the radius r h . 
resJts as shown in Fig. 2(d) v. ore obtained. One can readily sec that the ..ieasured data for the actual living 
T-re «< which is shown in Ftg 2(a) is in good agreement with the caJcutated data *x me phantom of the 
present invention which is shown in Fig. 2(d). and the agreement is particularly good at scattering angles in 
C-e range of 0 * 90 * . 

,io ideally, the phantom of the present invention is composed of polystyrene particles having many 
different values of radius. In practice, however, it is difficult and economically infeasible to obtain 
co-istyrene particles having many different values of radius and this necessitates a certain kind of 
aacrcximation. Thus, noting the two peaks indicated by A and B in Fig. 2(b). one may form a phantom by 
rrrcrg polystyrene particles having two different radii. Stated more specifically, polystyrene particles having 

is a radius r* ( = 0.025^0.001 urn) are mixed with polystyrene particles having a radius r 2 ( = 0-5 10. 01 5 urn), 
w-en these two types of polystyrene particles were suspended in water at respective concentrations of Ci 
(- 6.7%) and C2 ( = 2.5%). a phantom having the light-scattering characteristics snown in Fig. 3 was 
caained. The mean scattering cosine, g. of this phantom was 0.839. which is in very good agreement *ith 
r-ts value (g - 0.9) for the actual living tissues. 

20 The mean scattering cosine, g. is a parameter for the directivity of the scattering of particles and is 
eo-essed iry: 
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wrere e is the scattering angle and P(fl) is the scattering intensity. In the case cf complete forward 
scattering cne scattering angle 9 is 0* for all particles), g » 1: in the case of complete isotropic scattering. 
5 = 0: and in the case of complete backwaid scattering (d • 180* for alt particles), g - -l. 
30 A modification of the embodiment described above is explained below. While the distribution of 
scattering intensity was discussed in the foregoing embodiment, there exists another parameter to 
osscribed the light scattering, i.e.. coefficient of scattering us. The coefficient of scattering us is a 
parameter expressed by the reciprocal of the me-n optical pathlength which light travels b*^veen 
successive scatterings, and is given by: 
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•rare r fc is the radius of a particle. C k is the concentration of particles having the radius r k . and S k is the 
scattering cross section (as determined by Mie theory of light scattering) of the particles of the rac" - r k . 
Trarefore. d the relative cone nation of particles is determined as shown in Fig. 2(b) by performing 
operations on the matrix equation (1). the absolute concentration of the particles that yield a predetermined 
*s caefficieni of scattering u, can be determined. 

Besides use in evaluating the performance of a certain diagnostic device, the phantom of the present 
iruention has many other applications and three of them are briefly described be tow. 

so <v Determining optical path lengths in living tissues: 

When one wants to measure the change in the quantity of oxygen in the brain, he has to determine the 
concentrations of hemoglobin, oxidized hemoglobin and other vectors of oxygen transport by measunng the 
crage « fcght absorption at various wavelengths. In the absence of light scattering, the change in fight 
55 ataorption. A/ -ill be expressed by equation (4) to be given below. When, as shown in Rg. 4. incident 
Mgtt l« is taunc.<*d into the brain 1 . producing transmitted Bgnt I. AA is given by: 
a* - t • - • A C - Alog (Wo) (4) 

rare l is the length of optical path in the brain, a is the absorption coefficient of hemoglobin, and AC is 
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the change in oxygen concentration. Since I and U are measurable values and o is known, the change in 
oxygen concentration AC can be determined by equation (4). However, in ire presence of tight scattering, 
the effective value of optical pathiength is increased (several times the vate shown in Rg. 4). making it 
impossible to calculate AC by equation (4). If. in this case, an experiment is conducted using the phantom 
s of me present invention which has similar optical characteristics to the living assues. the effective value of I 
can be determined, thereby allowing AC to be calculated by equation (4). 



(2> Evaluation of imaging sucn as oxygen distribution: 

TO 

in imaging a certain phenomenon such as oxygen distribution in living assues. experimentation wttn a 
phantom and the evaluation of its results are indispensable. For this puraose. three light absorbers are 
contained in the phantom of the present invention as shown in Rg. 5. when* the phantom is indicated by 2 
and the light absorbers by 3. and the phantom 2 is illuminated with inoaert Bght By investigating the 
ts light transmitted through the phantom, the distribution of the intensity of traramrtted lig-nt car. be obtained in 
association with the light cfesoroers. and this allows exact evaluation of d o agnostic device of interest and 
precise diagnosis experiment, including, e.g.. the evaluation of image resofcflcn. 



7C (3) Analysis of light distribution in living tissues: 

In performing POT (photodynamic therapy), it is essential to know are extent of light propagation 
through living tissues of interest. This information can be obtained in an easy and exact way by using the 
phantom of the present invention. 

» Various modifications are possible with the present invention. For instarce. physiological saline may be 
used instead of water as a medium in which polystyrene particles are to be suspended. If a tight-absorbing 
dye is added in a suitable amount, the absorption characteristics oi are phantom can be adjusted 
independently of its scattering characteristics. 

As described in detail on the foregoing pages, a phantom having sirmar optical characteristics to living 

x tissues of interest can be produced in the present invention by merely incarcerating different-sized particles 
m a medium. Particularly high stability and reproducibility can be teaksec by selecting particles whose 
specific gravity is generally the same as that of the medium (for instarce. by suspending polystyrene 
particles m water). 



Claims 



1. A phantom having similar optical characteristics to living tissues with respect to light-scattering 
characteristics. compns.ng: 
«c a selected medium (2); and. . . 

a mixture of particles (3) suspended in said medium, of different-valued radus r k (k is a positive integer oi l 
to m) at respective values of concentration C k . wherein 

said values of concentration C t are determined on the basis of the following rnatnx equation. 
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where P^(fl.) »s a scattering intensity with said living tissues at a scattering angle 8. ft is a positive integer 
of 1 to n) and P(6..r*> is a scattenng intensity with said particles havrfic, said radius r* at said scattering 

^ 2. A phantom according to claim 1. wherein said values of concentraon C* are determined further on 
the basis of the following equation: 
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PS = t °k Sk 

5 where u« is a redproc* . iiean optical pathionpti and S» is a scattering cross-section of said particles of 

radius t k . § 
3. A phantom according to claim 1 or 2. wherein said particles(3) have a specific gravity which is 

generally the same as that of said medium (2). 

TO 
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FIG. 1 
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FIG. 2(c) 
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FIG. 5 
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